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Measurements of the sphericity of primary charged particles in minimum bias proton-proton colli- 
sions at y/s = 0.9, 2.76 and 7 TeV with the ALICE detector at the LHC are presented. The observable 
is linearized to be collinear safe and is measured in the plane perpendicular to the beam direction us- 
ing primary charged tracks with pj > 0.5 GeV/c in |rj \ < 0.8. The mean sphericity as a function of 
the charged particle multiplicity at mid-rapidity (N c h) is reported for events with different pj scales 
("soft" and "hard") defined by the transverse momentum of the leading particle. In addition, the mean 
charged particle transverse momentum versus multiplicity is presented for the different event classes, 
and the sphericity distributions in bins of multiplicity are presented. The data are compared with 
calculations of standard Monte Carlo event generators. The transverse sphericity is found to grow 
with multiplicity at all collision energies, with a steeper rise at low N c h, whereas the event generators 
show the opposite tendency. The combined study of the sphericity and the mean pj with multiplicity 
indicates that most of the tested event generators produce events with higher multiplicity by generat- 
ing more back-to-back jets resulting in decreased sphericity (and isotropy). The PYTHIA6 generator 
with tune PERUGIA-201 1 exhibits a noticeable improvement in describing the data, compared to the 
other tested generators. 
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1 Introduction 

Minimum bias proton-proton collisions present an interesting, and theoretically, challenging subject for 
detailed studies. Their understanding is important for the interpretation of measurements of heavy-ion 
collisions, and in the search for signatures of new physics at the Large Hadron Collider (LHC) and Fer- 
milab. However, the wealth of experimental information is currently poorly understood by theoretical 
models or Monte Carlo (MC) event generators, which are unable to explain with one set of param- 
eters all the measured observables. Examples of measured observables which are not presently well 
described theoretically include the reported multiplicity distribution [1-3], the transverse momentum 
distribution [4] and the variation of the transverse momentum with multiplicity [5-7]. 

In this paper, we present measurements of the transverse sphericity for pp minimum bias events over a 
wide multiplicity range at several energies using the ALICE detectors. Transverse sphericity is a momen- 
tum space variable, commonly classified as an event shape observable [8]. Event shape analyses, well 
known from lepton collisions [9-11], also offer interesting possibilities in hadronic collisions, such as the 
study of hadronization effects, underlying event characterization and comparison of pQCD computations 
with measurements in high Ej jet events [12-14]. 

The goal of this analysis is to understand the interplay between the event shape, the charged particles 
multiplicity, and their transverse momentum distribution; hence, the present paper is focused on the 
following aspects: 

- The evolution of the mean transverse sphericity with multiplicity for different subsets of events 
defined by the transverse momentum of the leading particle; 

- the behavior of the mean transverse momentum as a function of multiplicity; 

- the normalized transverse sphericity distributions for various multiplicity ranges. 

The results of these analyses are compared with event generators and will serve for a better understanding 
of the underlying processes in proton-proton interactions at the LHC energies. 

2 Event shape analysis 

At hadron colliders, event shape analyses are restricted to the transverse plane in order to avoid the bias 
from the boost along the beam axis [12]. The transverse sphericity is defined in terms of the eigenvalues: 
X\ > of the transverse momentum matrix: 



Since Sxy is quadratic in particle momenta, this sphericity is a non-collinear safe quantity in pQCD. For 
instance, if a hard momentum along the x direction splits into two equal collinear momenta, then the 
sum Y,i Pn will t> e half that of the original momentum. To avoid this dependence on possible collinear 
splittings, the transverse momentum matrix is linearized as follows: 




where (pxnPy,-) are the projections of the transverse momentum of the particle i. 
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The transverse sphericity is defined as 

s _ 2X 2 ( 
X2 + X\ 

By construction, the limits of the variable are related to specific configurations in the transverse plane 



This definition is inherently multiplicity dependent, for instance, Sj — > for very low multiplicity events. 

3 Experimental conditions 

The relevant detectors used in the present analysis are the Time Projection Chamber (TPC) and the Inner 
Tracking System (ITS), which are located in the central barrel of ALICE inside a large solenoidal magnet 
providing a uniform 0.5 T field [15]. 

The ALICE TPC is a large cylindrical drift detector with a central membrane maintained at -100 kV 
and two readout planes at the end-caps composed of 72 multi-wire proportional chambers [17]. The 
active volume is limited to 85 < r < 247 cm and —250 < z < 250 cm in the radial and longitudinal 
directions, respectively. The material budget between the interaction point and the active volume of the 
TPC corresponds to 1 1% of a radiation length, averaged in | r\ \ < 0.8. The central membrane divides the 
nearly 90 m 3 active volume into two halves. The homogeneous drift field of 400 V/cm in the Ne-C0 2 -N 2 
(85.7%-9.5%-4.8%) gas mixture leads to a maximum drift time of 94 jUs. The typical gas gain is 10 4 [7]. 

The ITS is composed of high resolution silicon tracking detectors, arranged in six cylindrical layers 
at radial distances to the beam line from 3.9 to 43 cm. The two innermost layers are Silicon Pixel 
Detectors (SPD), covering the pseudorapidity ranges |tj| <2 and |tj| <1.4, respectively. A total of 9.8 
millions 50 x 425 jltm 2 pixels enable the reconstruction of the primary event vertex and the track impact 
parameters with high precision. The SPD was also included in the trigger scheme for data collection. The 
outer third and fourth layers are formed by Silicon Drift Detectors (SDD) with a total of 133k readout 
channels. The two outermost Silicon Strip Detector (SSD) layers consist of double-sided silicon micro- 
strip sensors with 95 jltm pitch, comprising a total of 2.6 million readout channels. The design spatial 
resolutions of the ITS sub-detectors (a r< p x a,) are: 12 x 100 jim 2 for SPD, 35 x 25 /1m 2 for SDD, and 
20 x 830 jltm 2 for SSD. The ITS has been aligned using reconstructed tracks from cosmic rays and from 
proton-proton collisions [16]. 

The VZERO detector consists of two forward scintillator hodoscopes. Each detector is segmented into 
32 scintillator counters which are arranged in four rings around the beam pipe. They are located at 
distances z = 3.3 m and z = —0.9 m from the nominal interaction point and cover the pseudorapidity 
ranges: 2.8 < r\ < 5.1 and —3.7 < r\ < — 1.7, respectively. The beam-related background was rejected 
at offline level using the VZERO time and by cutting on the correlation between the number of clusters 
and track segments in the SPD. 

The minimum bias (MB) trigger used in this analysis required a hit in one of the VZERO counters or in 
the SPD detector. In addition, a coincidence was required between the signals from two beam pickup 
counters, one on each side of the interaction region, indicating the presence of passing bunches [1]. 

4 Data analysis 

MB events at y/s = 0.9 and 7 TeV (recorded in 2010) and at y/s = 2.76 TeV (recorded in 2011) have 
been analyzed using about 40 million events, each at 7 and 2.76 TeV, and 3.6 million at 0.9 TeV. Since 
no energy dependence is found for the event shape observable, we present mostly results for 0.9 and 7 




"pencil-like" limit 

1 "isotropic" limit 



TeV. 
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The position of the interaction vertex is reconstructed by correlating hits in the two silicon-pixel layers. 
The vertex resolution depends on the track multiplicity, and is typically 0.1 —0.3 mm in the longitudinal 
(z) and 0.2 — 0.5 mm in the transverse direction. The event is accepted if its longitudinal vertex position 
(z v ) satisfies \z v — zo\ < 10 cm, where zo is the nominal position. 

To ensure a good resolution on the transverse sphericity, only events with more than two primary tracks 
in |t]| < 0.8 and pj > 0.5 GeV/c are selected. The cuts on v\ and pi ensure high charged particle track 
reconstruction efficiency for primary tracks [7]. These cuts reduce the available statistics to about 9.1, 
4.2 and 0.42 million of MB events for the 7 TeV, 2.76 TeV and 0.9 TeV data, respectively. 

At 7 TeV collision energy, the fractions of non-diffractive events after the cuts are 99.5% and 93.6% 
according to PYTHIA6 version 6.421 [18] (tune PERUGIA-0 [19]) and PHOJET version 1.12 [20], 
respectively. In the case of single-diffractive events the fractions are 0.3% and 4.8%, while the double- 
diffractive events represent 0.2% and 1.6% of the sample as predicted by PYTHIA6 and PHOJET, re- 
spectively. 



4.1 Track selection 

Charged particle tracks are selected in the pseudorapidity range |tj| < 0.8. In this range, tracks in the 
TPC can be reconstructed with minimal efficiency losses due to detector boundaries. Additional quality 
requirements are applied to ensure high tracking resolution and low contamination from secondary and 
fake tracks [7]. A track is accepted if it has at least 70 space points in the TPC, and the % 2 per space 
point used for the momentum fit is less than 4. Tracks are rejected as not associated to the primary 
vertex if their distance of closest approach to the reconstructed event vertex in the plane perpendicular 
to the beam axis, do, exceeds 0.245 + (pj in GeV/c, do in cm). This cut is tuned to select primary 

Pt 

charged particles with high efficiency and to minimize the contributions from weak decays, conversions 
and secondary hadronic interactions in the detector material. 



4.2 Selection of soft and hard events 

The analysis is presented for two categories of events defined by the maximum charged-particle trans- 
verse momentum for |tj| < 0.8 in each event. This method is often used in an attempt to characterize 
events by separating the different modes of production. It aims to divide the sample into two event 
classes: a) events dominantly without any hard scattering ("soft" events) and b) events dominantly with 
at least one hard scattering ("hard" events). Figure 1 shows the mean transverse sphericity versus maxi- 
mum pj of the event obtained from minimum bias simulations at y/s = 7 TeV using the particle and event 
cuts described previously. Note that PYTHIA6 simulations (tunes: ATLAS-CSC [21], PERUGIA-0 and 
PERUGIA-20 1 1 [22] ) exhibit a maximum around 1.5 — 2.0 GeV/c, while PHOJET shows an intermediate 
transition slope in p™ x = 1 — 3 GeV/c. This observation motivated the choice of the following separation 
cut: "soft" events are defined as events that do not have a track above 2 GeV/c, while "hard" events are 
all others. The aggregate of both classes is called "all". The selection of 2 GeV/c has been motivated 
in the past as an accepted limit between soft and hard processes [23]. For parton-parton interactions 
the differential cross section is divergent for pj — > 0, so that a lower cut-off is generally introduced in 
order to regularize the divergence. For example in PYTHIA6, the default cut-off is 2 GeV/c for 2 — > 2 
processes. 

Table 1 shows the ratio of "soft" to "hard" events for ALICE data and the generators: PHOJET, PYTHIA6 
(tunes ATLAS-CSC, PERUGIA-0 and PERUGIA-2011) and PYTHIA8 version 8.145 [24]. It illustrates 
the difficulties to reproduce the evolution of simple observables with collision energy. 
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p™* (GeV/c) 



Fig. 1: Mean transverse sphericity versus /?™ ax for MC simulations at y/s = 7 TeV. Results are shown for PHOJET 
and PYTHIA6 (tunes ATLAS-CSC, PERUGIA-0 and PERUGIA-201 1) simulations. The events are required to 
have more than 2 primary charged particles in |tj| < 0.8 and transverse momentum above 0.5 GeV/c. 





0.9 TeV 


2.76 TeV 


7 TeV 


ALICE (data) 


5.7 


3.54 


2.36 


PHOJET 


8.53 


4.34 


2.52 


ATLAS-CSC 


10.95 


5.76 


3.41 


PERUGIA-0 


5.6 


3.26 


2.06 


PERUGIA-2011 


6.78 


3.64 


2.29 


PYTHIA8 


7.28 


3.92 


2.37 



Table 1: Ratio of the number of "soft" to "hard" events for data and MC generators according to the event selection 
criteria defined in the text. Corrections for trigger and vertexing inefficiency have been applied resulting in < 2% 
systematic uncertainty. 

4.3 Corrections 

The MC simulations used to compute the correction include transport through the detector and full re- 
construction with the same algorithms as the data. 

To correct the measured mean sphericity for efficiency, acceptance, and other detector effects, and to ob- 
tain it as the number of charged particles (A^h) in |t] | < 0.8 two steps were followed. First, the measured 
sphericity distributions in bins of measured mid-rapidity charged particle multiplicity (7Y m ) are unfolded 
using the detector sphericity response matrices. The unfolding implements a % 2 minimization with reg- 
ularization [25]. Second, to account for the experimental resolution of the measured multiplicities, the 
mean values of the unfolded distributions ((5 , x) unf ) are weighted by the detector multiplicity response, 
R(N c h,N m ). This procedure can be seen as 

(St) (^ch) = £ (S T ) unf (N m )R(N ch ,N m ) . (2) 

m 

Figures 2 and 3 show an example of the sphericity response matrix with a measured multiplicity of 25 
charged particles at mid-rapidity and the multiplicity response matrix, respectively. The MC simulations 
are based on the PYTHIA6 tune ATLAS-CSC. Different simulations were tested, and all produce the 
same results to within 1%. 

The sphericity distributions in four bins of multiplicity: (a) 3 < < 10, (b) 10 < N c ^ < 20, (c) 20 < 
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Fig. 2: Example of the sphericity response matrix for a measured multiplicity of 25 charged particles at mid- 
rapidity. The events are generated using the PYTHIA6 tune ATLAS-CSC (pp collisions at y/s = 7 TeV) and then 
transported through the detector. Particles and tracks with |tj| < 0.8 and pj > 0.5 GeV/c are used. 
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Fig. 3: Example of the multiplicity response matrix. The events are generated using PYTHIA6 tune ATLAS-CSC 
(pp collisions at s/s = 7 TeV) and then transported through the detector. Particles and tracks with |tj| < 0.8 and 
pj > 0.5 GeV/c are used. 



N c h < 30 and (d) 7V ch > 30 are also presented. The normalized spectra give the probability of finding an 
event with certain sphericity at given multiplicity. The normalized spectra were corrected bin-by-bin as 
follows 



P(St) \ atNch =P(S^) \ atNm xCi x C 2 , (3) 

where P(Sj) \ a tN m is the measured probability of finding an event with sphericity Sj in a bin of measured 
multiplicity (N m ). This probability is corrected by C\ and C%, which are computed using MC. C\ is the 
correction of the spectra at the measured multiplicity bin 

C ~ ^ I (4^ 
~ P(S m ) m ' 

and C2 corrects the probability by the migration from high to low multiplicity 

„ P(Sj) \atN ch 

^l^ T J \atN m 
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Fig. 4: Performance of the procedure to correct the reconstructed mean pj as a function of multiplicity for "all" 
events. The method is tested using PHOJET as input and applying corrections derived from PYTHIA. The MC 
true (PHOJET result at generation level) is compared with the corrected result after simulation and reconstruction. 



Contribution 


All 


Soft 


Hard 


Track selection cuts 


0.3% 


0.3% 


0.3% 


Event generator dependence 


0.5% 


0.5% 


2% 


Different run conditions 


1.0% 


1.0% 


1.0% 


Secondary track rejection 


< 0.8% 


< 0.8% 


< 0.8% 


Pile-up events 


0.2% 


0.2% 


0.2% 


Method (N ch < 5) 


< 5.0% 


< 5.0% 


< 11.0% 


Method (N ch > 5) 


< 1.5% 


< 1.5% 


< 1.5% 


Detector misalignment 


negl. 


negl. 


negl. 


ITS efficiency 


negl. 


negl. 


negl. 


TPC efficiency 


negl. 


negl. 


negl. 


Beam-gas events 


negl. 


negl. 


negl. 


Total (N ch < 5) 
Total (No, > 5) 


< 6.0% 

< 2.2% 


< 6.0% 

< 2.2% 


< 12.0% 

< 3.0% 



Table 2: Contributions to the systematic uncertainties on the mean transverse sphericity (St)- 



In the expressions, P(Sj) is the probability of finding an event with true sphericity S l T , where "true" refers 
to the value obtained at generator level. S { T and 5 T nf are the true and unfolded sphericity distributions, 
respectively. The latter are the results of the unfolding of the simulated measurements, i.e. PYTHIA6 
(tune PERUGIA0) corrected by PHOJET and vice versa. 

Finally, to determine (pT)(N C h), we take the mean pj by counting all tracks that pass the cuts discussed 
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Multiplicity range 




3-9 


10-19 


20-29 


Method 


< 


0.1% 


< 2.0% 


< 5.0% 


Event generator dependence 


< 


5.0% 


< 1.0% 


< 1.0% 


Pile-up events 


< 


1.0% 


< 1.0% 


< 4.0% 


Total 


< 


5.1% 


< 2.4% 


< 6.5% 



Table 3: Systematic uncertainties on the sphericity distributions, 
above as a function of measured multiplicity (N m ). Once we get (pj) m (N m ) we follow the approximation 

(Pt) (N ch ) = £ (pj) m (N m )R(N ch ,N m ) . (6) 

m 

Note that in this case an unfolding of the mean pi is not implemented. Figure 4 illustrates the perfor- 
mance of the procedure using PHOJET simulations as input. The response matrices are computed as 
above using the PYTHIA6 event generator. The corrected points are compared with MC at generation 
level. The differences, at high multiplicity, reach about 1.5%. 

4.4 Systematic uncertainties 

The systematic uncertainties on (St) are evaluated as follows. To minimize the adverse effects of pile-up 
on the multiplicity only runs with a low probability of multiple collisions were used. The parameter used 
to measure the pile-up level is the median of the Poisson distribution which is based on the recorded beam 
luminosity, it is assumed to characterize the probability to have n interactions reconstructed as a single 
event. Furthermore, a cut on the number of extra vertices reconstructed by the SPD was introduced. The 
systematic uncertainty was estimated from the differences between results using runs with the smallest 
and largest pile-up probability, and for the "all" sample was found to be less than 0.2%. The uncertainty 
due to the rejection of secondaries was estimated by increasing their contribution up to ~ 8%. This is 

done by varying the cut on the distance of closest approach (do > 0.0350+ 0,0 ^ Q , pj in GeV/c, do in 

Pi 

cm) of the considered track to the primary vertex in the plane perpendicular to the beam. The event 
generator dependence was determined from a comparison of the results obtained when either PYTHIA6 
or PHOJET were used to compute the correction matrices, and found to be of the order of few %. The 
most significant contribution to the systematic uncertainties is due to the method of correction. It was 
estimated from MC by the ratio true-Sx to corrected-Sr as a function of multiplicity. For example, 
the largest uncertainty is at low multiplicity (A^h ~ 3) for the "hard" sample, where it reaches ~ 11%. 
Different sets of cuts were implemented in order to estimate the systematic uncertainty due to track 
selection. Table 2 summarizes the systematic uncertainties on (St). In addition, other checks were 
performed to ensure an accurate interpretation of the results. For instance, when applying the analysis to 
randomized events (where the track azimuthal angles are uniformly distributed between and 2n), we 
obtain results that are about 10% larger than in data. The conclusion is that measured sphericity in data 
is not the result of a random track combination. Also, the analysis was applied to events with sphericity 
axes in different regions of the TPC, to ensure that the results are not biased by any residual geometry 
effects. 

In the case of the mean transverse momentum as a function of multiplicity the systematic uncertainties 
are taken from [7], the only difference being the method of correction. The uncertainty was estimated by 
applying the correction algorithm to reconstructed events generated with PYTHIA6, while the correction 
matrices were computed using events generated with PHOJET The final distributions were compared 
with the results at generator level. For the "all" sample the uncertainty reaches 1.5%, while for "soft" 
and "hard" it reaches 1.0% and 5.1%, respectively. 

For the case of the sphericity distributions in intervals of multiplicity, the main uncertainties are listed in 
Table 3. They were estimated following similar procedures as described above. 
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Fig. 5: Mean transverse sphericity as a function of charged particle multiplicity. The ALICE data are compared 
with five models: PHOJET, PYTHIA6 (tunes: ATLAS-CSC, PERUGIA-0 and PERUGIA-2011) and PYTHIA8. 
Results at ^/s = 0.9 and 7 TeV are shown in the top and bottom rows, respectively. Different event classes are pre- 
sented: (left) "soft", (middle) "hard" and (right) "all" (see text for definitions). The statistical errors are displayed 
as error bars and the systematic uncertainties as the shaded area. The horizontal error bars indicate the bin widths. 
Symbols for data points and model predictions are presented in the legend. 



5 Results 

In this section the results of the analyses are presented along with predictions of different models: PHO- 
JET, PYTHIA6 version (tunes: ATLAS-CSC, PERUGIA-0 and PE RUGIA-201 1) and PYTHIA8. 

5.1 Mean sphericity 

The mean transverse sphericity as a function of N c h at ^/s = 0.9 and 7 TeV is shown in Fig. 5 for the 
different event classes. The mean sphericity (right panel) increases up to around 15 primary charged 
particles, however, for larger multiplicities the ALICE data exhibit an almost constant or slightly rising 
behavior. For "soft" events and yfs = 0.9 TeV, the models are in agreement with the ALICE measure- 
ments over the full range of multiplicity, except for PYTHIA8 prediction, which is 5 — 10% lower. There 
is insufficient statistics to perform the unfolding for N c h > 18. At 7 TeV, the differences between models 
and data are below 10% for "soft" events. For the "hard" events, PHOJET, ATLAS-CSC, PERUGIA-0 
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Fig. 6: Mean sphericity versus multiplicity for (left) "soft", (middle) "hard" and (right) "all" events for y/s = 0.9, 
2.76 and 7 TeV. The statistical errors are displayed as error bars and the systematic uncertainties as the shaded 
area. 



and PYTHIA8 predict a lower (St) than observed in data, actually the differences between models and 
data are larger than 10% for multiplicities below 10 and larger than 40, that is true at 0.9 and 7 TeV. The 
differences observed are larger than the systematic and statistical uncertainties. It is interesting to note 
that PERUGIA-2011 describes the data quite well. The fraction of "soft" and "hard" events in data and 
MC simulations as a function of N c h (integral values given in Table 1) is found to be different between 
data and the event generators. At large N c i„ the event generators generally produce more "hard" events 
than observed in data. This difference is reflected in the "all" event class, since more "hard" events con- 
tribute in the case of the generators, while more "soft" events in the case of data. The largest isotropy in 
the azimuth is found at high multiplicity, > 40 (|tj| < 0.8, pt > 0.5 GeV/c), in a similar multiplicity 
region where the CMS collaboration discovered the long-range near-side angular correlations [26]. Com- 
paring the results at 0.9 and 7 TeV, it is seen that except for Pythia8 the predictions of models describe 
better the 0.9 TeV data than the 7 TeV ones. Lastly, the mean sphericity evolution with multiplicity at 
the three measured energies are shown in Fig. 6 for "soft", "hard" and "all" events at y/s = 0.9, 2.76 and 
7 TeV. The functional form of the mean sphericity as a function of N c h is the same at all three energies in 
the overlapping multiplicity region. 



5.2 Mean transverse momentum 

The mean transverse momentum as a function of N c h at yfs = 0.9 and 7 TeV is shown in Fig. 7. As seen 
in left panel, PERUGIA-0, PERUGIA-2011 and PYTHIA8 are within the systematic uncertainty bands 
of the data for soft events, though PYTHIA8 has a different functional form than the data. For the "hard" 
events there is a significant difference between the data and the generators above a multiplicity of about 
20, in particular for the 7 TeV data. For lower multiplicities, ATLAS-CSC has an overall different shape 
than the other generators. For "all" events, at 0.9 TeV PERUGIA-0 and PERUGIA-201 1 best reproduces 
the data, while the rest of the models do not give a good description. At 7 TeV, the calculations exhibit 
a change in the slope around N c h = 30, which is not observed in the data. At similar multiplicities, the 
MC mean sphericity reaches a maximum before it decreases with increasing multiplicity (Fig. 5). The 
similarity in the multiplicity dependence between (St) and (pj) suggests that the models may generate 
more back-to-back correlated high pj particles (jets) than present in the data. 
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Fig. 7: Mean transverse momentum versus multiplicity. The ALICE data are compared with five models: PHO- 
JET, PYTHIA6 (tunes: ATLAS-CSC, PERUGIA-0 and PERUGIA-201 1) and PYTHIA8. Results at ^ = 0.9 and 
7 TeV are shown in the top and bottom rows, respectively. Different event classes are presented: (left) "soft", 
(middle) "hard" and (right) "all". The gray lines indicate the systematic uncertainty on data and the horizontal 
error bars indicate the bin widths. 



5.3 St spectra in multiplicity intervals 

To disentangle the ambiguities between pj, Sj and multiplicity, the normalized transverse sphericity 
spectra (the probability of having events of different transverse sphericity in a given multiplicity interval) 
are computed at 7 TeV for four different intervals of multiplicity: N c h = 3-9, 10-19, 20-29 and above 
30. These are shown in Fig. 8 along with their ratios to each MC calculation. In the first multiplicity bin 
(A^ch = 3-9), the agreement between data and MC is generally good, but in the second bin (7V c h = 10-19) 
the ratio data to MC is systematically lower for Sj < 0.4 except for PERUGIA-201 1. In the last bin of 
multiplicity the overproduction of back-to-back jets (in the azimuth) reaches a factor of 3, and there is 
an underestimation of isotropic events by a factor 2. As in previous cases, the best description is done by 
PERUGIA-201 1. 

To obtain information about the interplay between multiplicity and (pj) through the event shapes, we also 
investigated the (pj) as a function of (St) in intervals of multiplicity. The study is presented using MC 
generators at = 7 TeV, but the conclusion also holds at the other two energies. Figure 9 shows (pj) 
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Fig. 8: Sphericity distributions in four bins of multiplicity: (upper-left) 3 < A^ c h < 9, (upper-right) 10 < iV C h < 19, 
(bottom-left) 20 < A^ c h < 29 and (bottom-right) N c h > 30 at ^/s = 7 TeV. The statistical errors are displayed as 
error bars and the systematic uncertainties as the shaded area. 



as a function of Sj for two multiplicity bins (top panels) along with the contribution of each sphericity 
bin (bottom panels) to the final (pj), i. e. the (pj) weighted by the value P(Sj). There are two points 
to emphasize. First, a large dependence of (px) on sphericity is observed for high multiplicities while 
at low ones the dependence is weaker. Second, the sphericity distribution determines the mean pi in a 
specific bin of multiplicity. For instance, for St = 0.3-0.4 PHOJET and ATLAS-CSC have nearly the 
same value of (pj), while the contribution to (pj) in the multiplicity bin is twice larger for PHOJET 
compared to ATLAS-CSC. Hence, the reproduction of the sphericity should be taking into account in the 
tuning of the MC generators. 
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The ALICE Collaboration 




Fig. 9: Mean pj (top) as a function of sphericity for two multiplicity bins (left) 3 < N c h < 9 and (right) No, > 30 
for minimum bias pp collisions at s/s = 7 TeV simulated with four different MC generators: PHOJET, PYTHIA6 
(tune ATLAS-CSC and PERUGIA-201 1) and PYTHIA8. Also the contributions of the different event topologies 
to the averaged mean pj are presented (bottom). 



6 Conclusion 

A systematic characterization of the event shape in minimum bias proton-proton collisions at = 0.9, 
2.76 and 7 TeV is presented. Confronted with the persistent difficulties of event generators to reproduce 
simultaneously the charged particle transverse momentum and multiplicity, the transverse sphericity is 
used to provide insight into the particle production mechanisms. The observables are measured using 
primary charged tracks with pj > 0.5 GeV/c in |tj| < 0.8 and reported as a function of the charged 
particle multiplicity at mid-rapidity (-/V c h) for events with different scales ("soft" and "hard") defined by 
the transverse momentum of the leading particle. The data are compared with calculations of standard 
Monte Carlo event generators: PHOJET, PYTHIA6 (tunes: ATLAS-CSC, PERUGIA-0 and PERUGIA- 
201 1) and PYTHIA8 (default MB parameters). 

The MC generators exhibit a decrease of (St) at high multiplicity with a simultaneous steep rise of (pr). 
On the contrary, in ALICE data (St) stays approximately constant or slightly rising (Fig. 5) accompanied 
with a mild increase in (pj) (Fig. 7). The mean sphericity seems to primarily depend on the multiplicity 
and not on ^/s (Fig. 6). At high multiplicity (N c h > 30) the generators underestimate the production 
of isotropic events and overestimate the production of pencil-like events (Fig. 8). It seems that the 
generators tend to produce large multiplicity events by favoring the production of back-to-back high-£>T 
jets (low St) more so than in nature. The level of disagreement between data and generators is markedly 
different for "soft" and "hard" events, being much larger for the latter (Figs. 5-7). It is worthwhile to 
point out that PERUGIA-2011 describes the various aspects of the data generally quite well, except 
for the mean pj, which it overestimates at high multiplicities. Our studies suggest that the tuning of 
generators should include the sphericity as an additional reference. 
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